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The 2-point angular correlation function w{9) (2PACF), where 9 is the angular separation be¬ 
tween pairs of galaxies, provides the transversal Baryon Acoustic Oscillation (BAO) signal almost 
model-independently. In this paper we use 409,337 luminous red galaxies in the redshift range 
2 = [0.440, 0.555] obtained from the tenth data release of the Sloan Digital Sky Survey (SDSS 
DRIO) to estimate 9bao{z) from the 2PACF at six redshift shells. Since noise and systematics can 
hide the BAO signature in the w — 9 plane, we also discuss some criteria to localize the acous¬ 
tic bump. We identify two sources of model-dependence in the analysis, namely, the value of the 
acoustic scale from Cosmic Microwave Background (CMB) measurements and the correction in the 
0BAO (z) position due to projection effects. Constraints on the dark energy equation-of-state param¬ 
eter w{z) from the 9bao{z) diagram are derived, as well as from a joint analysis with current CMB 
measurements. We find that the standard ACDM model as well as some of its extensions are in 
good agreement with these 9bao{z) measurements. 


I. INTRODUCTION 


Baryon Acoustic Oscillations (BAO) arise due to the 
competing effects of radiation pressure and gravity in the 
early Universe W- At z ~ 1000 photons and baryons 
decouple and a characteristic scale rg, corresponding to 
the sound horizon at the drag epoch, is imprinted in the 
maps of galaxy distribution and in the power spectrum 
of cosmic microwave background anisotropies^. In re¬ 
cent analysis, a tiny excess of probability to find pairs of 
galaxies separated by a distance equal to this comoving 
acoustic radius was revealed in the 2-point correlation 
function (2PCF) of galaxy catalogs, where it appears as 
a bump. The first detections of this BAO signature were 
obtained from galaxy clustering analysis of the Two De¬ 
gree Field Galaxy Survey (2dFGRS) [B] and from the 
Luminous Red Galaxies (LRGs) data of the Sloan Dig¬ 
ital Sky Survey (SDSS) [7j. A subsequent joint analy¬ 
sis of 2dFGRS and SDSS data yielded a BAO distance 
measurement with aggregate precision of 2.7% at redshift 
z = 0.275 [8] (see also i)- More recently, higher-z mea¬ 
surements at percent-level precision were also obtained 
using deeper and larger galaxy surveys 0 [Mil] (see 
also [131 for a recent review). 

The BAO signature defines a robust standard ruler, 
providing independent measures of the angular diam¬ 
eter distance Da{z) and the Hubble parameter H{z) 
through the transversal and radial BAO modes, respec- 
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1070 is the redshift of the drag epoch and Cs{z) is the sound 
speed of the photon-baryon fluid. 


tively. However, the detection of the BAO signal through 
the 2PCF makes necessary the introduction of a fiducial 
cosmology to calculate the radial distance to the cosmic 
objects then the comoving distance between them. On 
the other hand, the calculation of the 2-point angular cor¬ 
relation function (2PAGF) involves only the angular sep¬ 
aration between pairs, yielding model-independent infor¬ 
mation about 0 bao(-z) or, equivalently, Da{z), provided 
that robust estimates of the comoving acoustic scale are 
obtained. Therefore, determining the angular position of 
the BAO bump at several narrow redshift shells allows us 
to test the observational viability of different dark energy 
models through the 0bao ~ diagram. In this type of 
analysis, narrow redshift shell are necessary in order to 
avoid contributions from the radial BAO signal. 

In recent years, efforts have been made to map increas¬ 
ingly large volumes of the sky and to explore the cosmo¬ 
logical consequences of the BAO signature imprinted in 
galaxies distribution. One of the major multi-filter imag¬ 
ing and spectroscopic redshift surveys is the Sloan Digital 
Sky Survey (SDSS) that has been operating over fourteen 
years El- In this work we use the tenth public data re¬ 
lease (DRIO) [T5| of BOSS experiment, part of the SDSS- 
III project, to perform almost model-independent cos¬ 
mological analyses using the BAO transversal signature 
obtained from the 2PAGF. We use only 2-dimensional 
information, which restricts us to focus on the 2PACF 
in a set of narrow redshift shells. As the 2PAGF is usu¬ 
ally noisy due to systematic effects, we also introduce a 
double-tool methodology to identify the BAO bump. Our 
almost model-independent approach is possible due to 
some characteristics of the SDSS DRIO which were absent 
in previous releases. As will be clear in the next section, 
the most important one is the high galaxy number den¬ 
sity which allows it to have sufficient galaxy-correlated 
pairs revealing the angular BAO signature in small red- 
shifts shells {5z = 0.01 and 0.02). 

This paper is organised as follows. In Sec. H we de¬ 
scribe the basic methodology of our study, including the 
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relevant underlying equations for the 2PACF analyses. 
Sec. Ill discusses the observational data employed in this 
paper and the criteria for selecting the redshift shells. A 
new methodology to identify the BAO signal from the 
2PACF curves as well as the underlying cosmological as¬ 
sumptions to extract 0 bao(-z) from the LRGs catalogs 
are described in Sec. IV. Assuming a time-dependent 
parameterisation for the dark energy equation of state, 
in Sec. V we discuss the cosmological constraints on the 
parameters of dark energy models from the 0 bao('Z) di¬ 
agram. Our overall conclusions are summarised in Sec. 
VI. 


II. THE 2-POINT ANGULAR CORRELATION 
FUNCTION 


The two-point correlation function (2PCF), ^(s), s be¬ 
ing the comoving separation, is usually used to determine 
the BAO feature given a set of cosmic tracers like, e.g.. 
Luminous Red Galaxies (LRGs). This function is defined 
as the excess probability of finding two pairs of galaxies 
at a given distance and is obtained by comparing the 
real catalog to random catalogs that follow the geometry 
of the survey. Among several estimators of the 2PGF 
discussed in the literature (see, e.g., [TBUSO] !. one of the 
most commonly used is the one proposed in Ref. [20] (and 
adopted in our analyses): 

DD{s) — 2DR{s) + RR{s) 

=- mi) -^ 


Here, DD{s) and RR{s) correspond to the number of 
galaxy pairs with separation s in real-real and random- 
random catalogs, respectively, whereas DR{s) stands for 
the number of pairs with comoving separation s calcu¬ 
lated between a real-galaxy and a random-galaxy. As 
mentioned earlier, the comoving distance s between pairs 
of galaxies is calculated assuming a fiducial cosmology. 
In fact, assuming a flat universe, as indicated by recent 
GMB data |2T1122| , the comoving distance s between a 
pair of galaxies at redshifts Zi and Z 2 is given by 

s = \/r2(zi) -h r^{z 2 ) - 2 r{zi)r{z 2 )cos 0 i 2 , (2) 

where 0 i 2 is the angular distance between such a pair 
of galaxies, and the radial distance between the observer 
and a galaxy at redshift Zi, r(zi), depends on the cos¬ 
mological model adopted. For instance, for a flat dark 
energy model with equation-of-state parameter w(z), it 
reads 


c dz 

Hq Jq •y/Hm(l -|- z)^ -|- (1 — f2m)^(z) 


(3) 


where 


l{z) = exp 


3 f [I-t-w(0')dln (1-I-z')] 
. Jo 


Hq and flm are, respectively, the Hubble constant and 
the present-day matter density parameter. 

Analogously to the 2PCF, the 2PACF is defined as 
the excess joint probability that two point sources are 
found in two solid angle elements dHi and dH 2 with 
angular separation 0 compared to a homogeneous Pois¬ 
son distribution [23]. This function can be used model- 
independently, considering only angular separations in 
narrow redshift shells of small 6z in order to avoid con¬ 
tributions from the BAO mode along the line of sight 
(radial signal). The 2PAGF is now calculated as a func¬ 
tion of the angular separation 0 between pairs at a given 
redshift shell, i.e.. 


^^e).£mzJEm±3m, , 4 ) 

where the transversal signal of the acoustic BAO scale 
manifests itself as a bump at certain angular scale Ofit- 
Naturally, once the position of such BAO bump is lo¬ 
calised, it is possible to build the 0bao{z) diagram, i.e., 
the evolution of the angular diameter-distance with red¬ 
shift (see Sec. VI). 

One can obtain the expected 2PAGF, we, in terms 
of the expected 2PGF, ^e, considering a distribution of 
objects between redshifts zi and Z 2 as a function of the 
angular distance 9 between pairs PH US] , namely 


pOO pOO 

we{6,z)= / dzi4>{zi) / dz2 4>{z2)^E{s,z), 
JO JO 


(5) 


where z = (zi -|-Z 2 )/ 2 , with Z 2 = zi +5z, and 4>{zi) is the 
normalised galaxy selection function at redshift Zi. How¬ 
ever, if bin shells are narrow, 5z « 0, then zi « Z 2 and 
^e{s,Zi) ~ ^e{s,Z 2 )- Therefore, one can safely consider 
that ^e{s,z) depends only on the constant parameter z, 
instead of on the variable z (for shells of arbitrary 5z, 
see [26]). The function ^_e(s, z) is given by 


^e(s,z) 



k^Jo(ks) 5^ Pm{k,z ), 


( 6 ) 


where jo is the zeroth order Bessel function, Pm{k,z) 
is the matter power spectrum and the bias factor has 
been set at 6 = I since it does not affect the BAO peak 
position (for a broader discussion on the plausibility of 
the bias assumption and non-linear corrections, we refer 
the reader to m and references therein). 


III. THE DATA SET 

The phase HI of the Sloan Digital Sky Survey [l4] . 
SDSS-HI, had several goals which include, among oth¬ 
ers, a better understanding of the mechanism behind cos¬ 
mic acceleration from the analysis of the BAO feature in 
the galaxy clustering. Four experiments have been op¬ 
erating to provide the DRIO m data sets: the Sloan 
Extension for Galactic Understanding and Exploration 
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redshift intervals 

number of LRGs 

Z 

6z 

0.440 - 0.460 

21,862 

0.45 

0.02 

0.465 - 0.475 

17,536 

0.47 

0.01 

0.480 - 0.500 

40,957 

0.49 

0.02 

0.505 - 0.515 

21,046 

0.51 

0.01 

0.525 - 0.535 

22,147 

0.53 

0.01 

0.545 - 0.555 

21,048 

0.55 

0.01 


TABLE I: The six bin-redshift intervals and their properties; 
number of galaxies, mean redshift of the sample, z, and bin- 
width, Sz. Notice that contignous intervals are separated by 
a redshift interval of size 0.005 to avoid correlation between 
neighbours. 


(SEGUE-2), the Baryon Oscillation Spectroscopic Sur¬ 
vey (BOSS), Apache Point Observatory Galactic Evolu¬ 
tion (APOGEE) and Multi-Object APO Radial Velocity 
Exoplanet Large-area Survey (MARVELS). 

The BOSS experiment, as part of the SDSS-III project, 
mapped about 1.5 x 10® luminous red galaxies (LRG), as 
faint as i = 19.9, over 10,000 deg^, up to redshifts 2 < 0.7. 
BOSS used upgraded instruments to extend the redshift 
survey and map deeper than the SDSS-I and 11. This new 
spectroscopic experiment allows to map galaxy density 
with approximately n = 0.0002 — 0.0003 h®/Mpc^. 

The SDSS DRIO contains 409,337 LRGs in the north 
galactic hemisphere, with redshifts from z = 0.43 to 
z = 0.7, including new BOSS spectra. In order to de¬ 
tect transversal signatures as a function of redshift, we 
divided the data into six shells of redshift, as seen in Ta- 
bleH The selection of the width of the redshift shell is not 
an easy task: too narrow redshift shells may not contain 
enough galaxy correlated pairs to reveal the BAO signal, 
whereas if the redshift shell is too wide it decreases the 
angular BAO amplitude, mixing the transversal signa¬ 
ture with the radial contribution. After several analyses, 
we selected shells at: z = 0.45,0.47,0.49,0.51,0.53, and 
0.55, which are separated by a redshift interval of 0.005 
to avoid correlations between them. Details of these red¬ 
shift shells are presented in Table |l] 


IV. DETECTING THE BAO SIGNAL 

In previous works (see, e.g., [Ml HZ]), the 2PAGF 
was applied to galaxy surveys with photometric redshift 
(photo-z) data, where the redshift error depends on the 
range of wavelengths given by the filters. Unfortunately, 
the large error of the photometric redshift implies large 
uncertainties in the separations of the galaxies pairs. In 
addition to this problem, some analyses need a fiducial 
cosmology to derive cosmological parameters from the 
2PACF |M]- This strategy is usually adopted because it 
is very common to observe more than one single bump in 
any analysis of the 2PCF or 2PAGF, due to systematic 
effects present in the sample. Thus, when the 2PAGF 
shows more than one bump, one needs guiding principles 
to recognise the true acoustic scale. Since in our case 


we want to perform an analysis as model-independent as 
possible, we identify the angular-BAO signal only when 
a bump remains after applying the following double-tool 
methodology: 

1. Bin size criterium - Angular separations between 
pairs of galaxies are counted in bin intervals with 
bin width A0. Our analyses are performed for var¬ 
ious values of A6* in order to confirm if the BAO 
bump candidate persists or disappears. 

2. Small shifts criterium - We perform the 2PACF 
analysis changing the galaxies angular coordinates 
by small and random amount. Thus, the curve is 
smoothed and the bumps produced by systematic 
effects is removed. 

Our strategy is the following: if the BAO bump is 
present and robust, then it will remain after applying 
these criteria. If this happens, one can safely consider 
that the bump corresponds to a transversal BAO signa¬ 
ture. After that, the BAO bump is localised through 
a best-fit procedure, obtaining the OpiT value, which is 
then corrected to the 9bao value (corresponding to the 
ideal case in which 5z = 0) using a shift correction func¬ 
tion discussed in Sec. V. 


1. Bin size criterium 

In any histogram, like the 2PAGF, the choice of the 
bin size is a compromise between a noisy curve (when 
narrow bins are used), where possible signatures are hid¬ 
den by statistical noise and systematics, and a smooth 
curve (when wider bins are used), where possible signa¬ 
tures spread out and are almost invisible. Our leading 
criterion to decide the presence of a robust bump - like 
a BAO bump - is that it should persist in the 2PACF 
even when one changes the bin width, while noise bumps 
shall smear out or simply disappear. 

To optimize the choice of the bin width in our 2PACF 
analyses we consider several possibilities. Among them, 
we find three interesting cases for the angular interval, 
i.e., 9 G [0°,12°] (although the 2PAGF curves show 
only the interval of interest: [1.75°, 6.5°]). In Fig. [^we 
show the 2PAGF curves for the six redshift shells, con¬ 
sidering the cases where the number of angular bins is 
Af, = 31, 37, and 40. In each case the bin size is given 
by A6» = 12°/Aft. 

As an example, one can see in Fig. j^that for z = 0.55, 
the first bump located in d G [3°,4°] is present for the 
analysis considering 37 and 40 bins but its amplitude is 
significantly reduced (^ 30 %) when 31 bins are consid¬ 
ered. For comparison, the amplitude of the second bump, 
9 G [4°,5°], decreases only by ^ 6 % in these analyses, 
which shows its robustness relative to the first one, a clear 
indication that it corresponds to the acoustic scale. Sim¬ 
ilar results were also found for the other redshift shells. 
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FIG. 1: Each 2PACF plot assumes different redshift shells with almost 20,000 galaxies (see details in tableThe different 
9 bin corresponds to the interval from 0 to 12 degrees divided by Ni, = 31, 37, and 40 bins, i.e., /S.6 = 12.0/iVi,. To a better 
visualization of the 2PACF curves they were artificially shifted upwards. The error bars corresponds to one standard deviation 
obtained from the 50 random catalogs used in the analysis. 


2. Small shifts of the galaxies angular coordinates 

Various peaks or bumps can be observed in the 2PACF 
curves shown in Fig.[^ They contain information regard¬ 
ing not only the true BAO signature, but also to system¬ 
atic effects present in the distinct redshift samples used 
in our analysis. Distinguishing the BAO bump from sys¬ 
tematic bumps is the aim of the algorithm presented bel¬ 
low. The procedure is based on the hypothesis that the 
primordial BAO signature is present in the catalog. Ac¬ 
cording to this, we argue that if such signature is present, 
then it manifests itself as a robust bump at a given an¬ 
gular scale 9b AO- Instead, systematic effects, like signals 
produced by groups or clusters of galaxies, contribute 
with bumps at several angular scales, and do not survive 
to small perturbations of the galaxy positions. Conse¬ 
quently, one expects they will make no contribution to 
the 2PACF. 

Our algorithm follows two steps. First, we derive the 
2PACF corresponding to the case where no BAO signa¬ 
ture is present at all. This function is obtained averag¬ 
ing one hundred 2PACFs, each one obtained by chang¬ 
ing the angular positions of the galaxies by a random 
amount. The random displacements were performed fol¬ 


lowing Gaussian distributions with cr = 0.25,0.5, and 
1.0, which are equivalent to maximum displacements of 
9 ^ 1.25°, 2.5°, and 5°, respectively. In principle, this 
procedure can destroy any bump, resulting in a smooth 
averaged curve, as observed in the panels of Fig. (al¬ 
though similar results were obtained with N^, = 31 and 
37 bins, we used 40 bins in the analyses shown in Figs. 
and|^. In fact, one notices that for the cases with large 
values of cr, e.g., a = 0.5 and 1, all the peaks disap¬ 
pear. On the other hand, for cr = 0.25 only one peak 
remains, showing the robustness of the feature expected 
for the BAO signature. The second step is to compare the 
original 2PACF with the smoothed no-BAO curve corre¬ 
sponding to the case a = 0.5, case in which the noise is 
fully removed and comparison with the original 2PACF 
allows to identify the position of the BAO signal (the 
bin size criterium is also adopted to confirm this result). 
We then identify the BAO signature as corresponding 
to the excess in the original 2PACF with respect to the 
smoothed no-BAO curve, features that are more easily 
seen, e.g., in the shells z = 0.47 and z = 0.55. We 
emphasize that both criteria are applied to all redshift 
shells. 










5 



1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 


Q e 


FIG. 2: The continuous, dashed, and dot-dashed lines correspond to the 2PACF after changing the angular position of galaxies 
by a random amount following Gaussian distributions with a — 0.25,0.5, and 1.0, respectively, whereas the dots stand for the 
original data. In these plots we used Ni, = 40. 


A. Obtaining the 9fit values 

After finding the real BAO signature, we obtain the 
angular BAO scale using the method of Ref. [27] , which 
parameterises the 2PACF as a sum of a power law, de¬ 
scribing the continuum, and a Gaussian peak, which de¬ 
scribes the BAO bump, i.e., 

wfit [9) = A + B0'' + Ce , (7) 

where A,B,C,v, and apir are free parameters, 9f:t 
defines the position of the acoustic scale and a fit gives 
a measure of the width of the bump. If Sz = 0, the true 
BAO scale 9bao and 9fit would coincide. However, for 
6z ^ 0 this is no longer true because of projection effects 
due to the width of the redshift shells. Therefore, the 
angular correlation function given by Eq. ([^ has to be 
calculated for both 6z = 0 and 6z ^ 0 so that one can 
compare the position of the peak in the two cases. This 
will allows one to find a correction factor a that, given 
the value of 9fit found using relation Q, will provide 
the value for 9bao- 

In order to calculate the 2PCF given by Eq. Q, one 
needs the theoretical mater power spectrum, Pm(k,z). 


Models 



Wo 

Wa 

Ro“ 

Reference 

0.0226 

0.112 

-1 

0 

70 

Varying Q.ch'^ 

0.0226 

0.100 

-1 

0 

70 


0.0226 

0.140 

-1 

0 

70 

Varying state 

0.0226 

0.112 

-2 

0 

70 

equation 

0.0226 

0.112 

-0.8 

0 

70 


0.0226 

0.112 

-1 

1 

70 


0.0226 

0.112 

-1 

-1 

70 

Varying Ho 

0.0226 

0.112 

-1 

0 

65 


0.0226 

0.112 

-1 

0 

68 


0.0226 

0.112 

-1 

0 

72 


0.0226 

0.112 

-1 

0 

75 


“in units of km/s/Mpc 


TABLE II: Summary of the cosmological models parameters 
considered in our analysis. 


We use the CAMB software (Code for Anisotropies in 
the Microwave Background) jH]^, and assume a varying 


^ http://www.camb.info 
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FIG. 3: The 2PACF for six bin redshift intervals using the DRIO-SDSS data (bullets) and Eq. (continuous line). The 
amplitude of the BAG bump corresponds to C, the BAO location and the width are related to 9fit and a, respectively. In 
these plots we used Ni, = 40. 


« interval 

0 

a (%) 

Gfit (°) 

euz) n 

dsAO (°) 

O-BAO 

0.440-0.460 

0.45 

2.0815 

4.67 

4.96 

4.77 

0.17 

0.465-0.475 

0.47 

0.5367 

4.99 

4.77 

5.02 

0.25 

0.480-0.500 

0.49 

2.0197 

4.89 

4.60 

4.99 

0.21 

0.505-0.515 

0.51 

0.5002 

4.79 

4.44 

4.81 

0.17 

0.525-0.535 

0.53 

0.4847 

4.27 

4.29 

4.29 

0.30 

0.545-0.555 

0.55 

0.4789 

4.23 

4.16 

4.25 

0.25 


TABLE III: Estimates of 6bao{z) from SDSS DRIO LRG data. 


dark energy model with w(a) = wq + Wa(I — a), where 
a is the cosmological scale factor. In the present anal¬ 
ysis we assume a minimal model using six cosmological 
parameters 

{uJb,uJc,Q,T,As,ns}, ( 8 ) 

where uib = 0,bh? and Wc = ^ch'^ are, respectively, the 
baryon and cold dark matter densities, 0 is the ratio be¬ 
tween the sound horizon and the angular diameter dis¬ 
tance at decoupling, r is the optical depth to reionization, 
As is the overall normalization of the primordial power 
spectrum, and rig is the effective tilt. We consider purely 
adiabatic initial conditions, impose flatness and set neu¬ 
trino masses equal to 3.046 el^. We set the parameter 


values of the reference cosmological model used in this 
work as follows: ojbh'^ = 0.0226, ujch"^ = 0.112, 100 0 = 
1.04, T = 0.09, AsC^ = 2.2, Us — 0.96 with Hq = lOO/i 
km/s/Mpc. 

For each value z = the power spectrum P(k,z) 
was calculated and tabulated for values of k in the range 
[10“"^ — 2.38] h Mpc“^. From this table, we numerically 
calculate the integral in Eq. to obtain the spatial 
correlation function In order to minimize nu¬ 

merical uncertainties, we integrate analytically between 
neighboring points of the table by drawing straight lines. 
The fact that, in performing the numerical integration, 
we can not extend the upper limit of the integral to in¬ 
finity, gives rise to small oscillations that we smoothed 
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out using a standard moving average filter. The integral 
which gives the angular correlation function (Eq. was 
calculated using a top-hat distribution for the redshift 
shell selection function (f). We then numerically integrate 
we{S) for the redshift shells z and widths 5z listed in 
Table I and found the position of the peak 9^. By sim¬ 
ilarly proceeding and putting 6z = 0, we obtain 9^ and 
calculate the shift factor, defined as 


(0% - 

0 % 


( 9 ) 


Thus, the 9bao values will be given by the expression 


9bao{z,Sz) = 9fit{z) + a{z,Sz)9B{z) . (10) 


We calculated the shift factor a for several cosmologies 
(see Table II) in order to study its dependence upon the 
cosmological paremeters. Our overall conclusion is that 
the shift factor a is almost independent of the cosmolog¬ 
ical parameters in the range of values here considered, 
with the difference between 9fit and 9bao being < 2%. 
We also noticed that the small value of a is mainly due 
to the tiny redshift shell Sz < 0.02, chosen in our analy¬ 
ses, in that the lesser the value of Sz, the lesser the shift 
factor a is (for comparison, we refer the reader to Fig. 3 
of [57] ). The results for the reference model are shown in 
Table III. 


V. COSMOLOGICAL CONSTRAINTS 

In this Section we present cosmological parameter fits 
to the BAO data displayed in Table III. The angular scale 
9bao is related to the angular diameter distance Da{z) 
through 


dsAoiz) 


rs 

(1 -I- z)Da{z) ’ 


( 11 ) 


where Da{z) = r(z)/(l -I- z) and r{z) is given by Eq. 
([^. The evolution of 9bao with redshift, z, is presented 
in Figure 4 (left panel). The lines and coloured bands 
correspond to the standard ACDM cosmology, assuming 
ilm = 0.27, for two values of the comoving acoustic scale, 
i.e., Tg = 106.61 ± 3.47 /i“^Mpc and Ts = 100.29 ± 2.26 
/i“^Mpc given by the WMAP9 [5T] and Planck data [55] . 
respectively^. Along with the six data points obtained 
in our analysis, we also show other estimates obtained 
by translating the three-dimensional averaged distance 
parameter D^{z) = [(1 -I- z)'^D\{z)cz/H{z)Y^^ into an 
angular scale, with the fiducial cosmology used in each 


® An almost model-independent measurement of the BAO ruler 
length was recently reported in Ref. m- Using type la supernova 
and galaxy clustering data, the authors obtained Vs = 101.9 di 
1.9h-iMpc. 


reference. Such estimates are shown only for compari¬ 
son and will not be used in the statistical analysis that 
follows. 


As mentioned earlier, we consider a varying dark en¬ 
ergy model whose equation-of-state parameter evolves 
linearly with the scale factor, w(a) = wq -I- Wa(I — a) 
(w(a)CDM) and some of its particular cases, i.e., wq ^ 
— 1 and Wa = 0 (wCDM) and wq = — I and w^ = 0 
(ACDM). Plots of the resulting cosmological constraints 
are shown in Figure 4 (central and right panels). We 
find a tension between the 9bao{z) data derived in this 
analysis and the value of the acoustic scale Ts given 
by the Planck collaboration (see left panel of Fig. 4). 
For instance, assuming the wCDM cosmology, the com¬ 
bination of both data provides ilm = 0.39 ± 0.03 and 
Wo = -0.68 ±0.075. 


In what follows, we only display the results obtained 
assuming the acoustic scale estimate from WMAP9. 
Clearly, the 9bao data alone (blue contours) are con¬ 
sistent with a wide range of wq and w^ values, with the 
best-fit values being compatible with phantom scenar¬ 
ios in which w < —I. On the other hand, the combi¬ 
nation with the CMB data sharply limits the allowed 
range of w, favouring values of w ~ —1.0. This can be 
seen when we combine the BAO data points with mea¬ 
surements of the shift parameter (red contours), defined 
as TZ = yOm Hq/H{ z)dz, where zis is the redshift 

of the last scattering surface. To be consistent with 
the Ts value used in the BAO analysis, we use TZ = 
1.728 ± 0.016, also determined with model parameters 
set by WMAP9 data |5T|. The joint results (grey con¬ 
tours) improve significantly the cosmological constraints, 
providing 0^ = 0.29 ± 0.02 and wq = —0.91 ± 0.08 
(wCDM) and wq = —0.88 pto 0.20 and w^ = —0.22 ± 0.9 
(w(z)CDM) at 68.3% C.L. Assuming w = —1 (ACDM), 
we find 0^ = 0.334 ± 0.054 also at 68.3% C.L. 


Finally, it is important to mention the role of the 
acoustic scale in the 9bao analysis [see Eq. ([TI|)]. As 
discussed above, most of the estimates of this quantity 
currently available were obtained from CMB data (an 
exception being the measurement reported in Ref. [30] 1. 
From Eq. 0’ however, we can directly estimate 
from the 9bao data displayed in Table III, assuming 
a given cosmology. For the ACDM scenario, we find 
Ts = 101.2 ± 11.8 /i“^Mpc, which is in good agreement 
with both the WMAP9 and Planck values as well as with 
value obtained in [30j . In Fig. 5 we show the Dm - 
plane obtained when one assumes = 0.295 ± 0.034, 
as derived from current type la supernova data [JT]. In 
this case, we find Ts = 107.6 ± 2.3 /i“^Mpc, which is in 
full agreement with the ACDM estimate of Vg provided 
by the WMAP9 analysis and ~ 1.5a off from the value 
of the Planck collaboration. 





0BAO{°) 



FIG. 4: Left) The angular BAO scale as a function of redshift. As indicated in the figure, the blue data points correspond to 
the six measurements obtained in this paper (Table III) whereas the curves stand for the ACDM prediction with the acoustic 
scale fixed at the WMAP9 and Planck values. Central) Conhdence contours in the Llm - wq plane. Note that the combination 
between 6 'bao( 2) and CMB sharply limits the allowed range of the cosmological parameters. Right) The same as in the previous 
panel for the wq - Wa plane. 



FIG. 5: The flm - r^ plane obtained from the 9{z) data dis¬ 
played in Table III assuming the estimate of the matter den¬ 
sity parameter from current type la supernova data m- 


VI. CONCLUSIONS 

The baryon acoustic oscillations signal imprinted in 
the galaxy distribution is a key prediction of cosmologi¬ 
cal models, depending on the sound speed and expansion 
rate during decoupling. A decade after the first BAO 
detections, measurements of the BAO scale has become 
one of the main tools of precision cosmology which can 
be used to place sharp constraints on the main cosmolog¬ 
ical parameters, using the data from future large-volume 
galaxy surveys like, e.g., JPAS [33]. 

We have analysed the 2PACF of luminous red galaxies 
from the SDSS-DRIO data and measured the BAO signal 
in the distribution of galaxies at six redshift shells in the 
interval 0.440 < 2 < 0.555. Differently from the 2PCF 
analysis, the use of 2PACF involves only the angular sep¬ 


aration between pairs, yielding model-independent infor¬ 
mation about 0 bao('Z)- In practice, however, at least two 
relevant sources of model-dependence can be identified. 
First, due to projection effects {6z ^ 0), the true BAO 
scale Qbao differs from Ofit^ with the shift factor a be¬ 
tween them depending on the predicted power spectrum 
P(fc, z) of a given cosmology. For the range of cosmo¬ 
logical parameters displayed in Table II and the width 
of redshift shell considered in our analysis {5z < 0.02), 
we have found a ^ 2%. Nevertheless, as observed in 
Table m the largest shift obtained considering several 
cosmologies is 0.1°, which is one-third of the size of bin 
A9 in the 2PACF curves. In other words, this amounts to 
say that, in practical terms, the model-dependent shift is 
negligible. Second, in order to estimate the cosmological 
parameters from Eq. (Ill, an independent measurement 
of the acoustic scale is needed. However, as discussed 
in Sec. V, the current estimates of are not completely 
model-independent, with most of them being obtained by 
setting the ACDM parameters extracted from the CMB 
data. 


We have also introduced and applied a model- 
independent methodology to identify the true BAO 
bump. Such a procedure is different from what has usu¬ 
ally been done in the previous literature in which a given 
cosmology is taken as a guide to identify the true acoustic 
scale. After applying this methodology, we have derived 
a new dBAo('2^) sample, which increases signihcantly the 
current number of OBAoiz) data points available in the 
literature. From this 0bao{z) sample, we have derived 
cosmological constraints on a class of dark energy sce¬ 
narios with w(a) = Wq -b Wo(l — a). As shown in Fig. 
4 (central and right panels), although the 0 bao(-z) data 
alone are consistent with a wide range of wq and w^ val¬ 
ues, the combination with the CMB data sharply con¬ 
strains the allowed range of the dark energy equation of 
state, favoring values of w compatible with a cosmologi- 
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cal constant. 

A final aspect worth mentioning is the possibility of use 
the current ^baoC-z) sample to estimate some relevant 
cosmological parameters such as the present expansion 
rate, Hq, and the number and mass of neutrino species. 
In principle, this can be done by calculating the shift 
correction a, defined in Eq. (10), for a given cosmology. 
A detailed study exploring this possibility is currently in 
progress and will appear in a forthcoming communica¬ 
tion. 
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